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the Solutions of Its Nonlinear Differential Equations 
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A large-signal equivalent circuit for the Esaki or tunnel diode is presented that charac- 
terizes the dynamic and static response of this semiconductor device. Nonlinear differential 
equations are written on the model and a graphical-numerical solution technique is described. 
Analog computer solutions for the nonlinear equations are also discussed for various modes 
of operation. 



1. Introduction 

^^ost oF tJie equivalent circuit models presented 
lor the tuiuiel diode liave been small-signal types. 
The small-signal operation of this semiconductor 
device can be correctl}^ cliaracterized by its incre- 
mental resistance at a fixed bias point and some 
associated reactive elements. However, the large- 
signal departures from this bias point wliich typify 
bistable switching aiul relaxation oscillations make 
tJie model inadqeuate. In digital systems tlie active 
circuit elements a"e typicallx-^ reciuired to operate 
over extended ranges of their voltage-current charac- 
teristics. Tliis paper presents a large-signal equiva- 
lent circuit for the tunnel-diode (fig. 1). This 
circuit characterizes botli the static and dynamic 
responses of the device for all modes of operation. 
The network enclosed in the dashed box simulates 
tlie static V-I characteristic of the tunnel diode. 

Nonlinear differential equations are produced IVom 
the model and a graph ical-numei'ical solution tech- 
nique is described tluit is applicable to (hiven or 
self-excited systems. Tlie equations are readily 
solved by analog computer techniques, particularly 
when curve following is employed on a monotonic 
curve that is presented in this paper. 

An example of the graphical-numerical solution of 
the nonlinear equations for the i"elaxation oscillation 
case is presented in detail. Some comparisons are 
made between the graphical and aiudog computer 
solutions of bistable switching and sinusoidal oscilla- 
tion modes of operation. 

2. Tunnel Diode Large-Signal Equivalent 
Circuit and Solution Techniques 

Tlie validity of the small-signal tunnel diode model 
[1]^ has been extended over a large-signal range by 
incorporating the static terminal characteristics of 
the device into the model shown in figure 1. The 
lead inductance, Lg, the series body resistance, Rsj 
and the barrier capacitance, C, are assumed to m.ain- 
tain constant values over the operating range of 
interest [2]. Measurement techniques for these 
elements are described in the literature [3]. 



i Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Large-signm 
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tunnel diode equivalent circuit. 



This study was made using the GE ZJ-56 (1 ma 
peak ciuTeiit) tunnel diodes. Typical values for 
is, 11 si ^xnd Cm this device are: 

i, = 6.0 nh 

C^=5.0 pf 

Figure 2 sliow^s a static V-I plot of the terminal 
characteristics of the ZJ-56 tumuil diode. Curve A 
in figure 3 is a broken-lino approximation to this 
static plot which is simulated by the i{v) current 
generator and resistance, n, of figure 1 (dashed box). 
Tlie output from the i(v) current generator (see ap- 
pendix) is defined l)y either 



i(v)=ii{v)—ij{v) (fig. 3 curve B) 



(1) 



or by an analytic approximation to curve B 

i[v)=ki[l — exj) (-k2V^)] (fig. 3 curve C). (2) 

ii (v) is defined by the linear equation relating v and i 
in region 1 of curve A. ij{v) is defined by the linear 
equations relating v and i in each of the j linear regions 
(j= 1,2,3,4) of curve A. The resistance ri is set equal 
in value to the inverse slope of region 1 of curve A 
in figm^e 3. 

The following pair of equations can be written on 
the tunnel diode circuit in figure 4 with which all 
principal modes of operation such as switching and 
oscillating can be realized : 
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Figure 2. Static V-I characteristic for the GE ZJ-56 tunnel 
diode. 
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Figure 3. Curve A: Approximation to the GE ZJ-56 tunnel 
diode static V-I characteristic. 

Curve B: Broken-line approximation to the current from the 
i(v) generator. 

Curve C: Analytic approximation to the current from the 
i(v) generator. 
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cf^+^-kAl-exp i-k,v')]-H=0 (3) 
L^^+Ri,=E{t)-v. (4) 



E(t) is'the applied voltage, R=Rx+Rs and L=Li+Ls. 
Combining eqs (3) and (4) produces a nonlinear 
differential equation whose solution yields the junc- 
tion voltage, V, i.e., 



\dt 






--E{t)o>l 



(5) 




Figure 4. Tunnel diode circuit containing external Rj and L^ 
elements. 
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Figure 5. Typical construction increment in the F(v) versus 
V plane. 



where ooI=y7i' This equation does not lend itself 

to a closed form solution. It falls into the class of 
nonlinear equations of the type [4] 



dt' 



+m^^+g(v)=h(t). 



(6) 



After some investigation a graphical solution for this 
form was found and techniques for its application 
were worked out [5]. The solution process begins 
with an integration of eq (5) with respect to the time 
which yields 



f^+F(v) +£ g{v)dt=m) +K (7) 



where, 



^^"^ ^h £ [(^+¥)~^^'^^'' ^''P i-hv')'^ dv (8) 

9{v)=o^l [(l+f ) v-RkAl-^^^ (-^2»0]] (9) 

H{t)=o>l{' E{t)dt (10) 



i5L= integration constant. 
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Equation (7) is tlien converted into incremental 
form 

^+F{v;)+S„{t)^H{nM)~K (11) 



where, 



SAt)=^j[v{t)]dt. 



(12) 



A fixed value of M is chosen heuristically and an angle 

1 



a = arc tan 



A^ 



(13) 



is determined. It may be necessary to multiply 
equation (11) through by some constant in order to 
produce useful as. a is corrected to accommodate 
the departure of the abscissa-to-ordinate scaling of 
F{v) versus v from a one-to-one numerical relation- 
ship per unit length. Equation (1 1) is then rewritten 
into its final form preparatory to graphical solution as 

^v tan a+F{v^)=H{nM)+K-Sn{t). (14) 

The angle a is constructed ph^^sically on a template. 
Vertical lines of length H{nAt)—Sn{t)+KB.re con- 
structed in the F(v) versus v plane as the solution 
proceeds. The successive Av^s are determined by the 
construction technique shown in the typical segment 
in figure 5. Sn{t) is determined at each step of the 
process by applying the trapezoidal rule to the values 
of g{Vn) and g{Vn+i). This will be further detailed in 
the following section. 

3. Example of a Graphical-Numerical 
Solution of Equation (5), Relaxation 
Oscillation Case 

The circuit parameters were chosen to produce 
relaxation oscillations, i.e., 



R< 



L 



-r\ C 



and i?< |— r| 



where — r is the value of the inverse differential slope 
at each point in the negative resistance region of the 
V-I characteristic (|-r| ^15012). A type GE ZJ-56 
tunnel diode was employed in the circuit shown in 
figure 4. The values of the parameters were 

i=100Mh 
(7=10 pf 
E= 0.25«? (step) 
ri = 58.80 
]c,= 0.008 

co§=10^^ rad/sec. 

These immerical values were put into eq (5) which 
was then integrated in respect to time. In order to 
produce angle as of reasonable size the equation was 
multipUed through by 30X10"^ This yielded (in 
the form of equation (11)) 



Ax 



. ( ^^^^^^"' )+30[1.7t-.8[l-exp(-11.0^^)]] 

+30X10' r [1.34i'-.16[l-exp(-11.0z;2)]]rfi5 
=l.b{W){nAt)+K, (15) 



Since 



therefore 



AVn tan 



. /30><10-9\ 



a:=arctan 



/30X10"° \ 



(16) 
(17) 



For example, if A^ =30X10"^ sec is chosen as the 
construction increment then a =45°. The A^ chosen 
depends on the predicted response of the device. If 
At is too large, the construction will be extremely 
coarse. If At is too small, the solution will proceed 
very slowly. The value of A^ may be easily altered 
during tlie solution process in order to achieve 
optimum speed and resolution. The analytical error 
between the exact solution of eq (7) and its graphical 
solution is of the order of — l/2A^r^ in each interval. 
At this point the F[v) versus v and g(v) versus 
V plots were dra\\T:i where 



i^(?;) = 30[1.7^'-.8[l-exp(-11.0^'2)]] (18) 



and 



^(^0=30X10^ [1.340^-0.16 [l-exp(--11.0z;2)]] (19) 

Figures 6 and 7 A show these plots. 

Values of fl^(nAO = 7.5(10^) {nAt) n = l, 2, 3, . . . 
were tabulated. 

Since the F(^0 versus v curve was found to be 
reasonably plotted on a 10: 1 physical scale, the angle 
a was adjusted to compensate for this departure from 
a 1 : 1 relationship. This resulted m 

a=5°44' for A^ =30X10"^ sec 
a:=71°34' for A^^IXIO"^ sec 
a:=85°15' for A^=:0.3X10-' sec. 



Three templates were constructed containing these 
angles. 

The integration constant if, which is equal to v at 
^ = 0, was set equal to zero. The actual value of ^ 
was sufficiently small to make this completely 
reasonable. 

The Sn{t) term is evaluated at each step of the 
graphical solution by the trapezoidal rule. Values of v 
as a function of time such as v{nAt) and v[{n-\-l)At] 
are determined b}^ the successive constructions in the 
F{v) versus v plane. Then the values g[v{nAt)] and 
^b[(^+l)Ai^]] are directly obtained from the g{v) 
versus v plot. As the construction proceeds each 
quantity 
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Figure ^6. g(v) versus v plot employed in graphical solution of 
relaxation oscillator case. 



s„{t)- 



< 



giv[(nAt)]+g(v[(n+l)At]) \ 



At 



(20) 



is added on to the preceding total value of Sn{t) and 
this new total is subtracted from the next applied 
H(nAt). This gives the length of the vertical con- 
struction line, H(nAt)—Sn(t) (as shown in fig. 5). 

As each construction is made, the time required 
for each AVn to occur is known from the relationship 
of A^ to the angle a employed in the construction. 
When the H(nAt)—Sn(t) line falls below the F(v) 
curve the construction reverses direction and pro- 
ceeds towards the origin. 

The angle a was changed during the construction 
process in order to produce maximum resolution and 
solution speed. For example, in the jump region the 
angle a=85°15' (for A^-0.3X10"^ sec) was 
employed. 

Figure 7 A displays the H(jiAt) — Sn{t) locus in 
the F{v) versus v plane for the relaxation 
oscillation case. Figure 7B shows the junction volt- 
age, V, versus time that was derived from the graphi- 
cal construction. A single pulse from the actual 
ZJ-56 relaxation oscillator pulse train is shown super- 
imposed upon a pulse produced by the graphical 
method which was based upon the analytic approxi- 
mation of eq (2) . The graphical process will produce 
a train of similar pulses if it is continued. 
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Figure 7. (A) F(v) versus v plot employed in the graphical 
solution of the relaxation oscillator case {solid line) . 

The dashed line is the locus of n(nAt)—Sn(t) obtained during the graphical 
construction. 
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(B) Junction voltage, v, versus time obtained from graphical 
solution (solid line) superimposed upon actual pulse {dashed 
line) . 



4. Comparison of Analog and Graphical 
Results 

Equations (3) and (4) were meclianized on the 
National Bureau of Standards Mid-Century 500 
Analog Computer employing curve-following tech- 
niques for both the broken-line and analytic i(v)^s 
as defined by eqs (1) and (2). The method of con- 
structing the broken-line i(v) function is described 
in the appendix. Since these curves w^ere relatively 
smooth and monotonic, (as compared to curve A, 
fig. 3) the curve-following operated well. The R, 
L, and C elements were adjusted to produce bistable 
switching and (nearly) sinusoidal oscillations. 
Graphical solution technigues were then applied for 
these same operating conditions. Figure 8 displays 
the response of the system to a ramp driving function 
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Figure 8. (A) Graphical construction in the F(v) versus v plane. 

(B) Analog computer and graphical solution for junction v 
versus time {ramp-driven switching 7node). 
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Figure 9. (A) Graphical construction in the F(v) versus v 
plane. 

(B) Analog computer and graphical solution for junction v 
versus time {oscillatory mode). 
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FifiURp: 10. (A) Simulated bistable switching circuit with 
additional pulse triggering source. 



{B) Analog computer inechanization for circuit in figure 10 A. 



witli the circuit operating in the switching mode, 
i.e., J^>\-r\, /?-47012, L=5.0 nh and C=2.0 pf. 
Figure 9 shows the sinusoidal response to a voltage 
step with E<C\—r\,R=lG^, L=1 .b nhand C=2.0 pf. 
The bistable switching case was simulated on the 
analog computer for another circuit configuration 
(fig. lOA). Here the tunnel diode was steadily 
biased with a voltage, E, that produced a stable 
singuhir point near the peak in region 1 (curve A, 
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fig. 3). The circuit was then triggered over to the 
second stable singular point in region 4 by an injected 
current pulse, ip. It was possible to observe the 
effects of pulses of different duration and amplitudes 
on the triggering process. This technique can be 
used to study the optimum trigger requirements for 
the tunnel diode and other negative resistance 
devices. Figure lOB shows the analog computer 
mechanization. The i{v) function eq (1) was plotted 
with a conducting material and curve-followed 
during the simulation process. 

5. Appendix 

Crisson [6] has presented a method for producing 
a voltage-controlled negative resistance at the input 
terminals of an ''ideal amplifier'^ with feedback. It 
can be shown that the network enclosed by the 
dashed line in figure 1 displays the same character- 
istics as the * 'ideal amplifier.^' 

The "ideal amplifier" shown in figure 11 A can be 
represented by the circuit shown in figure llB. The 
resistance, Vj, seen looking into the 1-2 terminals of 
the circuit is 



rr 



R1R2 



R1+R2-M, 



(i= 1,2,3,4.) 



(21) 



where Mj is a variable mutual impedance term. 
When Af,>7?i+7?2 then rj<0 which produces a nega- 
tive resistance at terminals 1-2. Mj can therefore 
be the basis for the simulation of the tunnel diode 
static V-I characteristic with its alternating positive 
and negative regions. Resistance n (see fig. 1) is 
set equal to the inverse slope of region 1 in curve A 
figure 3. Therefore, if Mi is set equal to zero in 
this region, then the parallel combination of Ri and 
i?2 is chosen to be equal in value to Vi. In any one 
of the other three regions the variable Mj can be 
be determined from 



M,= (i2.+E,)(l-^)- 



(22) 



A conversion can be made from the circuit in 
figure llB to the circuit enclosed in the dashed box 
in figure 1 in which the current generator output is 



i(v)\j=ii 



<'>('-7) 



=^{v)-ijiv) 



(23) 



(24) 



where, ii(v) and ij(v) are defined by the functional 
relationship of the current to the voltage in each 
linear region of curve A in figure 3 and each region 
is bounded by the breakpoints on the curve. Curve 
B in figure 3 is the monotonic i(v) curve that is 
produced by eq (24). This curve is easily followed 
by curve-following analog techniques. An analytic 
approximation to curve B was also produced for use 
in obtaining a total differential equation for the 
junction voltage, v, of the tunnel diode. The 
approximation is 



(A) (e) 




(B) (D 



Figure 11. (A) Ideal amplifier. 
(B) Network equivalent to the ''ideal amplifier.'' 



i(v)=ki[l-exp (-k2v')] (fig, 3 curve C) (25) 

where, ^i=0.008 and k2=n.O for the GE ZJ-56 
tunnel diode type that is used in this study. 

The i{v) curve B for use in analog simulation is 
constructed by differentiating eq (24) in respect to 
V which gives 



di(v) 



dv 



__ dii(v) dij(v) 



dv 



dv 



(26) 



For each of the j linear regions (j =1,2,3,4:) of curve 
B a line with slope -^^ is drawn from region 



dv 
breakpoint-to-breakpoint starting at the origin. 

Note that ^ =0. 



dv 
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with the analog computer mechanization. 
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